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Here we thoroughly discuss some weak points of the thermal model 
which is traditionally used to describe the hadron multiplicities 
measured in the central nucleus-nucleus collisions. In particularly, 
the role of conservation laws, the values of hard core radii along 
with the effects of the Lorentz contraction of hadron eigen volumes 
and the hadronic surface tension are systematically studied. It is 
shown that for the adequate description of hadron multiplicities 
the conservation laws should be modified, whereas for the descrip- 
tion of hadron yield ratios the conservation laws are not necessary 
at all. Also here we analyzed the usual criteria for the chemical 
freeze-out and found that none of them is robust. A new chem- 
ical freeze-out criterion of constant entropy per hadron equals to 
7.18 is suggested and a novel effect of adiabatic chemical hadron 
production is discussed. Additionally, we found that the data for 
the center of mass energies above 10 GeV lead to the temperature 
of nil surface tension of about To = 147 ± 7 MeV. This is a very 
intriguing result since a very close estimate for such a tempera- 
ture was obtained recently within entirely different approach. We 
argue that these two independently obtained results evidence that 
the (tri)critical temperature of the QCD phase diagram is between 
140 and 154 MeV. 



1. Introduction 

Experimental data on heavy ion collisions has tradition- 
ally been described by the thermal model [1-27]. The 
thermal model core assumption is that fireball produced 
in the relativistic nuclear collision reaches thermody- 
namic equilibrium. Such an assumption allows one to 
describe the multiplicities of particles registered in the 
experiment using two parameters, namely the tempera- 
ture T and baryo-chemical potential The extracted 
values of T and fib not only decribe the experimental 
data, but they also give an essential information about 
the last stage of fireball evolution, when the inelastic col- 



lisions cease to exist, but the elastic collisions between 
hadron and the decay of resonances take place. This 
stage is usually called the chemical freeze-out. The ther- 
mal model was initially used for the AGS and SPS data 
[T3] and was subsequently employed to describe the data 
collected at SIS PUT5], SPS [TB] and RHIC [T7H2D]. Us- 
ing the thermal model it was possible to correctly predict 
the hadron ratios measured at LHC [2 , while the only 
wrong prediction for LHC was made for p/ir~ ratio |21| . 
An analysis of the energy dependence of thermal param- 
eters extracted from fits of the experimental data estab- 
lished the line of chemical freeze-out [22 . Consequently, 
thermal model is an established tool for particle produc- 
tion analysis and chemical freeze-out investigation. 

However, the thermal model suffers from several weak 
points which should be accounted for in more careful 
studies. The present paper is just devoted to a critical 
analysis of the thermal model and contains several direc- 
tions to develop and improve it . First of all we would like 
to notice that the term "thermal model" is a common 
name for a set of similar models, each having its spe- 
cific features such as the strangeness supression factor 
j s , the inhomogenious freeze-out scenario |14 | IT^ | l2^ | l?7] 
etc. Here we consider the minimal thermal model with 
two major parameters T and fif,, following the approach 
of Andronic et al., [2j, and below we briefly formulate 
the main problems of the model to be analyzed in the 
present work. 

• Particle table. In order to describe the exper- 
imental data using the thermal model one needs 
the masses, the widths and the decay branching 
ratios of all existing resonanses. In principle, the 
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mass spectrum of the hadronic resonances that are 
heavier than 2,3 GeV is known poorly and, hence, 
they could create a problem. However, recently 
it was shown that the large width of heavy reso- 
nances leads to their strong suppression [28] and, 
hence, their contribution into the thermodynamic 
functions of hadronic phase is negligible. 

Note that not only the parameters of hadrons from 
the "tail" of mass spectrum are poorly known. 
For example, both the mass and the width of 
cr(600) meson are not well established, but the 
thermal model predictions are strongly influenced 
by the values of the mass and width of this meson 
[3], while for many other baryons the width and 
branching ratios of decays are not well established 
at all. Thus, the particle table is one source of 
uncertainty of the thermal model. 

• Hard-core spheres radii value. The ideal gas 
description has proven to be unsatisfactory [4j long 
ago. The simplest way to introduce an interac- 
tion between hadrons is to use the repulsive hard 
core potential, since the attraction between them 
is usually accounted for via many sorts of hadrons 
|28| . In the simplest case this potential depends 
only on a single parameter - the hard core radius. 
In general case, each particle type may have its 
own hard core radius. But for the sake of conve- 
nience and simplicity the hard core radius is usu- 
ally taken equal for all particles. The usual value 
for such a radius is r — 0.3 fm. This value is 
motivated by the hard core volume known from 
nucleon-nucleon scattering [1 . There are, however, 
two restrictions on the hard core radii range: (i) 
they should be small enough to satisfy the condi- 
tion V e ig en <C V . i.e. that the total eigen volume 
of all particles V e i gen should be much smaller than 
total volume of the system V; (ii) on the other 
hand, these radii shouldn't be too small, because 
otherwise the model will lead to a contradiction 
with the lattice quantum chromodynamics (QCD) 
thermodynamics data [5]. Thus, there is a cer- 
tain freedom in defining the hard core radii values 
which, so far, was not systematically exploit to de- 
scribe the whole massive of existing experimental 
data. 

• Conservation laws. Also here we would like to 
discuss the baryon charge and isospin projection 
conservation laws. It was suggested to use them in 



the form [2] 

Si n iHi — hinit/V , 
J2i n i B i = Bimt/V . 

The initial values are chosen I^init = — 20 and 
Binit = 200 |2J, neglecting the fact that only the 
part of initial particles belong to the midrapidity 
region. Below we study the role of these conser- 
vation laws and show that such a treatment leads 
to physically unrealistic freeze-out volumes and to 
very bad description of hadron multiplicities, while 
the particle yield ratios (we use such a term for the 
ratio of multiplicities in order avoid a confusion) 
description can be extremely good. 

• Multiplicities fit. The fit of hadron multiplicities 
is usually performed using three parameters: T, 
and V [2j. Below we show that such a procedure 
combined with above mentioned conservation laws 
is mathematically ambiguous and it leads to the 
problems with the imposed baryonic charge con- 
servation law. 

In addition to these usual features of the thermal 
model we would like to thoroughly investigate the role of 
the Lorentz contraction of eigen volumes of hadrons to 
essentially improve the previous analysis [24 , 25 , 29J and 
to study the influence of the hadronic surface tension on 
the fit of the freeze-out parameters. As we argue the 
latter may provide us with new information about the 
critical temperature value of the QCD phase diagram. 
Such a comprehensive analysis of the different features 
of the thermal model and the experimental data will also 
allow us to elucidate the correct criterion of the chemical 
freeze-out which is a very hot topic nowadays. 

The work is organized as follows. The basic features 
of the thermal model are outlined in the next section. In 
Section 3 we discuss the important problems related to 
the conservation laws and propose their solutions. The 
discussion of the existing chemical freeze-out criteria and 
formulation of a new criterion along with a novel effect 
of the adiabatic chemical hadron production is given in 
Section 4, while Section 5 is devoted to the model refor- 
mulation for the multi-component hadron gas mixture. 
In Section 6 we investigate the values of the hard core 
radii and study the effect of the Lorentz contraction of 
hard core spheres. Section 7 is devoted to the analysis 
of the hadronic surface tension, while Section 8 contains 
our conclusions. 
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2. Model formulation 

In order to study the role of conservation laws in a form 
suggested in [2] and employed in their subsequent publi- 
cations it is, first of all, necessary to reproduce the results 
obtained in that work. For this purpose let us consider 
the Boltzmann gas consisting of s sorts of hadrons hav- 
ing the temperature T and the volume V . Each z-th sort 
is characterized by its own mass m, and chemical poten- 
tial /ii. Suppose that the number of particles of i-th sort 
is iVj. Then its canonical partition function is 



Z can (T,V,N u ...,N s ) = 

H- V I ... i ^ 



n 

i=l 



(2tt) 3 



exp 



T 



d A k 



N, 



(1) 



Here gi = 25 + 1 is the degeneracy factor of i-th hadron 
sort, k is the particle momentum. The corresponding 
grand canonical partition function reads as 



Z 



gr.can. 

X] ■■■ exp 

N 1= N h =Q 



T 



x Z can (T,N u ...,N s ). (2) 
From Q one gets the number of particles of each sort: 
Ni = 4>i{T,mi,gi)eTV 



9iV 
(2tt) 3 



exp 



T 



d A k. 



(3) 



Following the commonly accepted approach, we consider 
the conservation of baryon charge B, strangeness 5* and 
isospin projection J3 on average: 



TV 



^ riiSi = Sinn = , 
i=i 

J2 niBi = B mit /V = 200/V , 
i=i 

N 

Y, n ihi= himt/V = -2Q/V . 



N 



(4) 
(5) 
(G) 



These conservation laws define the value of the total 
chemical potential for the hadron of sort i as /i,; = 
Bi ■ fib + Si ■ fi s + ■ /i/ 3 , where the quantities Bi, 
Si and I$i denote, respectively, the baryonic, strange 
and isospin projection of such a hadron, while the corre- 
sponding chemical potentials are denoted as /if,, fj, s and 

Mis- 



The interaction of hadrons and resonances is usually 
accounted for by the hard core repulsion of the Van der 
Waals type [S] as 



P = Pi, 



.gas 



• exp 



p ■ b 



nl exp 



(-*) 



1 



pb 
T 



(7) 



where the pressure ptd.gas an d the i-th charge density 
rif" of an ideal gas is modified due to hard core repulsion. 
Here b — ^-(2R) 3 is the excluded volume for the hard 
core radius R, which in actual calculations was taken 
to be R = 0.3 fm for all hadrons. The usual Van der 
Waals correction affects the particle densities, but has no 
effect on particle ratios [S] . While its effect on the charge 
and particle densities may be strong, for the freeze-out 
densities obtained at and above the highest SPS energy 
\/Smn — 17.6 GeV the excluded volume correction leads 
to a reduction of the densities by about of 50 percent. 

The resonance decays is usually accounted for in the 
following way: the final multiplicity of hadron X consists 
of the thermal contribution N^ 1 and the decay ones: 



N 



Jin 



X 



N$ 



j\jdecay 



N$ 



^^fBr(Y 



N- 



where Br(Y — > X) is the decay branching of the Y-th 
hadron into the hadron X. The masses, the widths and 
the decay branchings were taken from the particle tables 
used by the thermodynamic code THERMITS |7|. 

The width T of the resonance of mean mass m is 
accounted for by replacing the Boltzmann distribution 
function in the particle pressure by its average over the 
Breit-Wigner mass distribution as 



exp 



-\/F 



poo 

J M (x~- 



T 



dxj 



d 6 k 



-m) 2 +r 2 /4 



/exp 



d 3 k 



dx 1 



(9) 



/■OO 

J M (xi-m) 2 +r 2 /4 

where M is the dominant decay channel mass. 
3. Role of conservation laws 

Using the thermal model formulated in previous section 
we fitted the hadron yield ratios in the energy range 
from AGS to RHIC, i.e. for VS/vjv = 2.7 200 GeV. 
We used the x 2 minimization for all the ratios available 
for this energy range as the fit criterion. The present 
consideration is very similar to that one used in |2 . The 
main sources of difference are listed below. 
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• The Boltzman statistics is used here instead 
of the quantum statistics employed in |2J. This 
allows us to essentially fasten the simulations since 
the momentum integration can be done only once 
for each hadron species. We checked that for the 
freeze-out temperatures T > 50 MeV obtained here 
the difference of the results due to the Boltzman 
statistics is almost negligible. 

• The charm conservation is not accounted for 
by the present model, since it is important only 
for the charmed particles multiplicities description 
which is not considered here. 

• The particle table used here is slightly different 
from that one of [2J, but this does not lead to a 
big difference in results. Although, in contrast to 
[2], we do not fit the mass and the width of c(600) 
meson. 

• Inclusion of the resonance width is done in 
this work for all values of colliding energy ^/Snn, 
while in |2J the width was accounted only for the 
AGS energy range. 

As it is seen from Fig. [T] the experimental hadron 
yield ratios are reproduced very well within the present 
model. The dependence of the chemical freeze-out fitting 
parameters on ^/Snn is given in Fig. U In Fig. § we 
show only the statistical errors for the obtained fit, while 
for the parameters of work [2] the shown errors account 
for the systematic and the statistical ones. As one can 
see from Fig. [2] the discrepancy between the results of 
the present model and that one of [2] is within the error 
bars. 

Similarly to [2 we found that both the chemical freeze- 
out temperature T and baryonic chemical potential fif, 
are almost independent of the initial value of the baryon 
charge Bi n a and the initial value of the isospin projec- 
tion Ia init . However, we found that the freeze-out vol- 
ume V which stands on the right part side of conserva- 
tion laws (|5| and ^ is very sensitive to them. The ob- 
tained chemical freeze-out volume dependence on >/Snn 
for B init = 200 and Iz init — —20 is shown in Fig. [3] 

From the larger symbols in Fig. [3] one can clearly see 
that for the center of mass collision energies y/S^N = 
2.7 - 4.3 and s/S^ = 12 - 200 GeV the found chem- 
ical freeze-out volume is so large that it exceeds the 
volume of kinetic freeze-out 0. Here we found that 
unlike the hadron yield ratios, the chemical freeze-out 
volume is very sensitive both to the excluded volume 
correction and to the values of parameters B init and 
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Fig. 1. The examples of the particle yield ratios description. The 
dots denotes the experimental values, while the lines show the 
result of fit. Upper panel: v / 5jviV = 6.3 GeV, T=139 MeV, fl b 
= 503 MeV, the mean square deviation per degree of freedom is 
X 2 /NDF = 4.8/4. Lower panel: v / 3iVJV = 130 GeV, T=169 MeV, 
H b = 31 MeV, X 2 /NDF = 3.4/9. 



Is inU . From Eq. ^ one can deduce that the larger 
value of the excluded volume b corresponds to the larger 
value of the chemical freeze-out volume V since larger b 
value one obtains the smaller particle concentrations nj 
and, consequently, the larger volume V = Bi n i t /^2 n i^i- 
Therefore, the minimal chemical freeze-out volume cor- 
responds to an ideal gas, i.e. for 6 = 0. This mini- 
mal chemical freeze-out volume is also shown in Fig. [3] 
Despite the absence of the excluded volume correction, 
chemical freeze-out volume for an ideal gas remains huge. 
Hence, we conclude that the initial values in conservation 
laws, namely Bi n i t and I^ init are of crucial importance 
for an extraction of chemical freeze-out volume. 

From the comparison of the hadron multiplicity N = 
n ■ V obtained by the thermal model and its experimcn- 
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Fig. 2. Dependence of thermal model fitting parameters on the 
center of mass collision energy v^Sjvjv- Upper panel: the chemical 
freeze-out temperature T vs. y/SpfN- Lower panel: the chemical 
freeze-out baryonic chemical potential /n, vs. \/Snn- The re- 
sults obtained from the fit of hadron yield ratios (circles) with the 
conservation laws and from the fit of hadron multiplicities (open 
squares) are compared with that ones obtained in [2] (crosses). 



tal value dN / dy\ y= o measured at zero rapidity, one can 
conclude whether the thermal model provides a reliable 
description of hadron multiplicities. Such a comparison 
for 7r + and K + mesons is shown in Fig. [4] Since the 
hadron yield ratios shown in Fig. 1 are described well 
by the thermal model, then one would expect that , if the 
multiplicity of a single hadron type is in a good agree- 
ment with the experiment data, then the multiplicities 
of all other sorts of hadrons should be well described too. 
However, from Fig. |4]one can see that at \/Snn > 10 
GeV the experimental values of multiplicities are much 
smaller than the theoretical ones. At lower %/ Snn there 
is no such a problem, despite the big chemical freeze-out 
volume V. Our conclusion is that for higher energies the 




VSnn, GeV 



Fig. 3. The chemical freeze-out volume vs. V-Siviv for the ideal 
hadron gas and the hadron gas with the hard core radii 0.3 fm. The 
smaller symbols correspond to the fit of hadron yield ratios with 
all the conservation laws j4|-{6]) accounted for, while the larger 
symbols are obtained by the fit of hadron multiplicities ignoring 
Eq. j5| (see text for details). 

value of parameters Bi n u and |/3 <n4t | should be taken 
smaller than for lower energies. In particular, according 
to Fig. [4] at -v/Sjvjv = 200 GeV the value of Bi n a and 
|/3 iriit | should be about 10 times smaller than tat ones 
for the low collision energies. 

A different way to describe the hadron multiplicities 
was used in [2 . The chemical freeze-out volume V was 
treated there as a free parameter. Let's show that such a 
treatement leads to mathematical ambiguity. Consider 
the conservation laws (|4])-(|6]) together with the following 
expression for system pressure 



p = p ld (T, fi b , fi s , fi h ) -exp 



p b 



(10) 



Such a system of equations has six unknowns, i.e. T, 
fibi A*s, (J-i 31 V, p, and four equations. Hence, two un- 
knowns should be treated as free fitting parameters. If , 
however, one treats three unknowns as free parameters, 
then one of the equations may be not satisfied in general. 
More specifically, if T, ^ and V are the free fitting pa- 
rameters, then the baryon charge conservation equation 
(j5| may be broken down. To demonstrate this explic- 
itly we have considered the thermal model fit with three 
free parameters - T, /if, and V and ignored the baryon 
charge conservation equation ([5]), while the isospin pro- 
jection conservation law (|6| was used to find the chemi- 
cal potential /i/ 3 . After fitting the experimental hadron 
multiplicities dN/dy\ y=0 (not their ratios!) for the same 
energy range as before, we found the resulting baryonic 
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Fig. 4. Upper panel: 7T" 1 " multiplicity at chemical freeze-out vs. 
V^NN- Lower panel: K + multiplicity at chemical freeze-out vs. 
VSnn- The circles correspond to the experimental data, whereas 
the stars are found from particle densities as N = n ■ V 

charge as Sb — V ■ Yli=i n i^i summing up the densities 
rii of all baryons and anti-baryons multiplied by their 
baryonic charge Bi. Clearly if Eq. (|5| is satisfied, then 
this sum Sb should match the value of Bi n n = 200. Fig. 
[5j however, demonstrates that Eq. (|5| cannot be satis- 
fied. Although the chemical freeze-out temperature and 
baryonic chemical potential obtained by the fitting of the 
hadron multiplicities do not differ essentially from that 
ones found by the fit of hadron yield ratios (see open 
squares in Fig. [2]), the freeze-out volumes obtained from 
the multiplicity fit are essentially smaller (see smaller 
symbols in Fig. [3| and more physically adequate for 
\/Snn > 5 GeV than that ones found from the fit of 
hadron yield ratios. 

Additionally, here we found that the values of Bi n i t 
and Ig (nit should strongly depend on collision energy: at 
VSnn = 200 GeV they are about ten times smaller than 




Fig. 5. The obtained baryonic charge Sb = V ■ ] fH^i (dots) 
vs. VSnn- The dots are calculated from the baryon charge 
conservation, whereas the line corresponds to the expected value 

Binit = 200. 

at the AGS energies. The hadron yield ratios are not 
sensitive to the values of Bi„i t and Is init , and hence the 
baryon charge and J 3 conservation can be neglected for 
the description of hadron yield ratios. If, however, one 
supposes that Bi n u = const > and Iz inU — const > 0, 
then the description of hadron multiplicities completely 
fails. Evidently, one can describe the hadron multi- 
plicities by introducing \/S~nn dependence of B init and 
Iz init values, if such dependences are known. Since such 
dependences are unknown, then one has to ignore the 
baryon charge ^ and isospin projection (JsJ) conserva- 
tion laws, and to fit the parameters T, /Zb, /j,j 3 and V to 
describe the hadron multiplicities or to fit the parame- 
ters T, a bl uj 3 in order to get the description of hadron 
yield ratios. Note, however, that the strangeness con- 
servation laws Q does not create such problems and, 
hence, it should be always obeyed. 



4. Chemical freeze-out criteria and adiabatic 
chemical hadron production 

The thermal model discussed above allows us to clarify 
the long standing question on the physically appropriate 
chemical freeze-out criterion which is widely discussed 
|22l |2"5] . The most popular chemical freeze-out crite- 
ria are (I) the constant value of the mean energyper 
hadron, (E)/(N) ~ 1.08 GeV, (II) the constant value 
of the entropy density to the cube of the temperature, 
s/T 3 ~ 7, and (III) the constant value of a total baryon 
and antibaryon density ng + n§ ~ 0.12 fm -3 . The cri- 
terion (I) is believed to be more robust, while the cri- 
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teria (II) and (III) show strong dependence on the hard 
core radius value [23]. We have performed the analysis 
and found the criteria (II) and (III) are not obeyed at 
all, while the criterion (I) validity depends essentially 
on the thermal model parameterization. The validity of 
these statements for the criteria (I) and (II) are, respec- 
tively, demonstrated in the middle and the upper panels 
of Fig. [6j Moreover, the thermal model results that we 
extracted from [2] are very similar to our ones, despite 
several differences in the parameterization of these two 
models. 

Our detailed analysis shows that there exists a much 
more robust chemical freeze-out criterion than all previ- 
ously discussed ones. This novel criterion corresponds to 
the constant value of the entropy per number of parti- 
cles which in terms of the entropy density s and hadron 
number density p par t can be expressed as 



Ppart 



7.18. 



(11) 



The lower panel of Fig. [6] shows that for two different 
parameterizations of the thermal model the ratio - s - 
stands between 6.6 and 7.6, i.e. the deviates withing 
± 8 % only, while the values of the center of mass energy 
of collision change for two orders of magnitude! Such a 
behavior of the — - — quantity evidences for the adia- 

Ppart 

batic chemical hadron production in heavy ion col- 
lisions. 



5. Multi-component gas and hard core radii 

Although it is clear that the hadron radii can serve as 
the parameters of the thermal model, they, however, are 
rarely treated as the free parameters. The common ap- 
proach is to fix one radius for all hadrons. The value of 
this radius was discussed [11J and it ranges from 0.2 fm to 
0.8 fm|12). However, the value r=0.3 fm that was taken 
from nucleon-nucleon scattering |10| seems to be an es- 
tablished value. To investigate the role of the hard core 
radii we introduce the different radii for mesons and for 
baryons R m and Rb, respectively. Also to study an influ- 
ence of the Lorentz contraction of hard core radii on the 
chemical freeze-out parameters T, fib and on the hadron 
yield particle ratios we need the hadron gas model which 
accounts for different values of their eigen volumes. For 
this purpose we use an approach developed in 1 2"oll2T)] . 
Below we give the necessary theoretical apparatus to 
study the multi-component hadron gas mixture, whereas 
the results of the global fit for the cases with and without 
Lorentz contraction are given in the subsequent section. 
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Fig. 6. Different chemical freeze-out criteria. Upper panel: s/T 3 
at chemical freeze-out vs. \ASjvjv- Middle panel: {E)/(N) at 
chemical chemical freeze-out vs. v'Sjvjv- Lower panel: the novel 
criterion of chemical freeze-out s/p ~ 7.18. The errors shwon for 
the results obtained here are combined the statistical and system- 
atic errors. The results of present work (squares) are very similar 
with that ones extracted from [2]. 
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Consider again the Boltzmann gas of s hardron species 
in a volume V at a temperature T. Let Ni be a quantity 
of the i-th sort of hadrons 



N = 



The total number of particles is M — J2i=i ^i- ^ is as- 
sumed that for every two sorts of hadrons i and j there 
is the excluded volume bij. Then one can introduce the 
excluded volume matrix B = (bij). Naturally, it is sup- 
posed that the matrix B is symmetric, i.e. b^ = bji. 

The canonical partition function can be obtained by 
adding the particles of some eigen volume one-by-one 
and taking in account all the corresponding excluded 
volumes of the previously added particles. Such an ap- 
proximation was suggested in |H] and it gives the fol- 
lowing expression for the canonical partition of the Van 
der Waals hadron gas mixture 



Zvdw(T, V, Ni) = 



n 



Ni 



NA 



V - 



N 1 BN 



-, M 



M 



(13) 



where the thermal particle density ^(T, m, g) is defined 
in p|, and N T is the transposed matrix to that one given 



by (121 



In the next step we write the grand canonical partition 
function (GCPF) as 



OO OO 



JVi=lW 2 =l N s = l \i=l 



T 



x Z 



VdW 



.(14) 



It is well known |3 1] that in the thermodynamic limit 
the GCPF can be replaced by the maximal term of the 
multiple sum in Z (the maximum term method). Sup- 
pose that the array N* gives the maximal term of Z. 
Then the system pressure is given by 



p/T = lim — 

1 " V^oo V 



lim — In 

V-s-oo V 



n 



.4 



N*\ 



x IV- 



(N*) T BN* 



M* 



,(15) 



where Ai — <j>i exp \_^f\ ■ Let us find N* from the maxi- 
mum conditions (i = l..s): 



d 



on* 



In 



n 



.4 



NT I 



V ■ 



(N*y BN* 



M* 



0. (16) 



Performing the differentiations, one gets 



& = Ai exp -£2^6j 



est 



(17) 



(12) with & 



V-- 



and 



(18) 



NT 



Using (JTtJ) one can express the hadron densities rii = ^ 
and the system pressure p as 



6 



J=i 



(19) 



The solution of the system of equations ( 17 )-( 19 1 defines 



the hadron densities for the multi-component gas. 
In a special case when all the elements of the excluded 



volumes matrix are equal 6y = vq Eqs. ( IT I — ( 19 1 give 



& = A l exp(-pu /r) , 

1 l+pvo/T ' 

^/r=(E:=i^)-cxp(-p« /T 1 ). 



(20) 



In this case the ratios of two particle densities from ( 20 1 



match that ones of the mixture of the corresponding 
ideal gases for an arbitrary value of vq, while the parti- 
cle densities themselves may essentially differ from the 
particle densities of the ideal gas. 

6. Results for hard core radii with the Lorentz 
contraction 

In this subsection it is assumed that all baryons have 
the same hard core radii Rb and all mesons have the 
same hard core radii R m . Then performing the global 
fit we would like, first, to find the pair of radii (R mi Rb) 
that provides the best fit, and, second, we would like to 
study the influence of Lorentz contraction of the chemi- 
cal freeze-out parameters. 

To simplify the numerics let us define that two hadrons 
belong to the same type, if their excluded volumes are 
equal. The number of equations in the system (17) is 
equal to the number of partcle types. Hence, the case 
with the Lorentz contraction included is more compli- 
cated, because instead of two sets of particles one should 
treat each hadron type as a new kind of particles. In or- 
der to avoid the large number of equations in system 
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(17), the particles heavier than 900 MeV are considered 



non-relativistically, i.e. they all belong to two sorts: ei- 
ther to the baryons with the hard core radius Rb or to 
heavy mesons with hard core radius R m . For two parti- 
cles i and j, both heavier than 900 MeV, the element of 
excluded volume matrix 6^ is ^(Ri + Rj) 3 - For other 
cases the second virial coefficients bij are calculated us- 
ing Eqs. (3) and (4) of [29 by the direct translation of 
one ellipsoid around the other with the subsequent aver- 
aging of the obtained excluded volume over the ellipsoid 
positions. 

The fit procedure is the same as described in the pre- 
ceding subsection: we fit the hadron yield ratios by T, 
Hb and V, respect the isospin projection conservation 
law Q to find the value of /i/ 3 , but ignore the baryonic 
charge conservation law (|5j. This our effort to study the 
role of the Lorentz contraction is inspired by the fact 
that the conventional thermal model has problems with 
the light mesons description. Thus, to describe the pion 
multiplicity and the ratios containing the pions it was 
proposed to introduce R n which is smaller than all the 
other hadron radii [5J[3n]. Another example of difficul- 
ties with the light mesons is a "strangeness horn", i.e. 
the peak in the K + /ir + ratio. Its description was finally 
improved by fitting the ct(600) meson mass and width 
[3], the obtained description is far from being very good. 

Let us see whether the Lorentz contraction might re- 
solve these problems. At a given temperature the eigen 
volume of the lighter particles decrease more than that 
one of the heavier ones. Consequently, the excluded vol- 
ume of lighter particles gets smaller compared to the ex- 
cluded volume of heavier ones |24 [[25ll2"9] . This provides 
us with the natural explanation of the fact that pion hard 
core radius R^ is smaller than other hard core radii. On 
the other hand, this might improve the strangeness horn 
description. It was also shown |29J that the Lorentz con- 
traction removes the causality paradox from the thermal 
model, i.e. at high densities the speed of sound does 
not exceed the speed of light. Therefore, it is necessary 
to incorporate the Lorentz contraction into conventional 
thermal model. 

To compare the models with and without the Lorentz 
contraction, we have chosen hard core radii Rb=0A5 fm, 
Rm = 0.3 fm and have found the new best-fit T and 
Hb values for the case with Lorentz contraction, see Fig. 
[7| From this fi gure one can conclude that the baryo- 
chemichal potential is almost unaffected by the Lorentz 
contraction, while the temperature is slightly higher for 
the case with the Lorentz contraction. 

It is interesting to check, whether the inclusion of the 
Lorentz contraction improves the strangeness horn de- 
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> 
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n no Lorentz contraction 
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O 
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Fig. 7. The chemical freeze-out parameters obtained from the fit 
with (circles) and without (rectangles) the Lorentz contraction. 
Upper panel: the chemical freeze-out temperature T vs. \/Sjvjv- 
Lower panel: the chemical freeze-out baryonic chemical potential 
Hb vs. V Sat jv- 



scription. From Fig. [8] we conclude that there is a small 
improvement, which is not enough to qualitatively im- 
prove the strangeness horn description. 

The important result is that the Lorentz contraction 
inclusion provides us with the better fit quality for any 
pair of radii (R m ,Rb). From Fig. [9] one can see the dif- 
ference between the x 2 /NDF values found without the 
Lorentz contraction and without it. Obviously, when 
both radii are small, then the correction due to the 
Lorentz contraction is small too. At (R m ,Rb) = (0.3, 
0.4) fm Ax 2 /NDF w 0.1, while X 2 /NDF itself is 1.48. 
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Fig. 8. Strangeness horn description improvement. 
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The simplest way to obtain the best fit radii is to per- 
form a global fit. including the radii into the fitting pro- 
cedure. However, in the (i? m ,i?f,) plane there exist the 
domains, where x 2 stays almost unchanged. For exam- 
ple, in the case without the Lorentz contraction \ 2 is the 
same along the line Rf, = R m , which follows directly from 
(20). This makes the straightforward global fit rather 
difficult. Therefore, we perform the fit procedure of par- 
ticle ratios in central nucleus- nucleus collisions at y/SwN 
= 2.7, 3.3, 3.8, 4.3, 4.9, 6.3, 7.6, 8.8, 12, 17, 130, 200 GeV 
for each pair of the radii {R m , Rb) and find the domains, 
where \ 2 differs from its minimal value less than 10%. 
The results are shown in Fig. |10| 



7. Determination of hadronic surface tension 

Recently the extremely important role of the surface ten- 
sion of quark gluon bags was realized within the exactly 
solvable models for the deconfinement phase transition 
with the tricritical |321f33j and the critical [34_ endpoints. 
It was shown [32-34 that the (tri)critical endpoint ap- 
pears due to vanishing surface tension coefficient, while 
at low baryonic densities the deconfinement phase tran- 
sition degenerates into a cross-over just due to the neg- 
ative values of surface tension coefficient. The existence 
of negative values of the surface tension coefficient at 
the cross-over temperature for demonstrated analytically 



10 



ISSN 2071-0186. Ukr. J. Phys. XXXX. Vol. XX, No. X 



INVESTIGATION OF HADRON MULTIPLICITIES AND HADRON YIELD RATIOS IN HEAVY ION COLLISIONS 



Collision energies set, VSjvjv 


X 2 /NDF without surface tension 


X 2 /NDF with surface tension 


ao.MeVfm- 2 


T , MeV 


2.7 - 7.6 


25.8135/33 = 0.782 


25.8043/31 = 0.832 


0.91 ■ 10" 2 


61 


2.7 - 200 


103.096/82 = 1.2573 


103.036/80 = 1.288 


-1.37- 10" 2 


57 


2.7 - 62.4 (no 130 and 200) 


85.51/65 = 1.3156 


85.268/63 = 1.3534 


-3.21 • 10~ 2 


62 


12, 17, 62.4, 130, 200 


62.5452/37 = 1.69 


62.1454/35 = 1.776 


0.654 


147 



Table 1. Results of the global fit, including the extracted surface tension parameters. 



within the model of color confining tube j35j . Using this 
model it was possible to predict the value of (tri)critical 



temperature of QCD phase diagram T a 



152.9 ±4.5 



MeV [36 using the plausible assumption on the tem- 
perature dependence of the surface tension coefficient 
u{T) — a 2 (\ — ^ [37] which is typical for ordinary 
liquids. 

Since the lattice QCD is not reliable at the non-zero 
values of the baryonic chemical potential it would be 
interesting to study the surface tension for hadrons at 
the chemical freeze-out. The surface free energy can be 
written as F sur f = a(T)S, where the hadron surface S 
is given by its hard core radius R as S — 4irR 2 . Inclu- 
sion of the surface free energy into the thermal model is 
equivalent to adding the term a(T)S to the total chem- 
ical potential. If S is the same for all hadrons, such 
a surface tension correction does not affect the hadron 
yield ratios. Therefore, we have taken i?f, = 0.45 fm, R m 
= 0.3 fm, to have the noticeable radii difference. In the 
actual simulations the surface tension coefficient cr(T) 
was parameterized as 

a (T) =a (l-~Y (21) 

Here oo and To > are the free parameters to be found 
from a global fit. Note that for T < T such a tem- 
perature dependence coincides with the famous Fisher 
droplet model parameterization [37 , whereas for T > T 
it is in line with the recent findings |32H36| . We have 
performed several global fits with the parameterization 
(21), using different data sets. The results are listed in 



the Table Q] 

From this table one can conclude that the inclusion 



of the surface tension in the form (21 ) does not improve 
the fit quality, however, since in all cases the value of 
X 2 /NDF is almost the same as without accounting for 
the surface tension, it also does not spoil the fit quality. 
This fact allows us to take the obtained values of pa- 
rameters (To and To rather seriously. It is not surprising, 
that the value of <7q is close to zero, otherwise the siz- 
able surface tension of hadrons could be already found. 
The really surprising fact is that for the center of mass 
energies \p5jss~N > 12 GeV the parameter T = 147 ± 7 



MeV is extremely close to the critical temperature value 
T ce p — 152.9 ± 4.5 MeV found in [35] more than a year 
ago using entirely different approach. Of course, the 
reason of why the global fit that includes the low energy 
data gives essentially lower value of the parameter T 
should be understood, and, hence, the investigation of 
the hadronic surface tension should be continued using 
both the experimental data on hadron production and 
the lattice QCD data. 

8. Conclusions 

In this work we performed a comprehensive analysis of 
the experimental hadron multiplicities within the ther- 
mal model. The considered thermal model has two hard 
core radii (Ri, for baryons and R m for mesons) and two 
new elements: inclusion of the Lorentz contraction of 
eigen volumes of hadrons and treatment of hadronic sur- 
face tension. Using this model we studied the role of the 
imposed conservation laws Q and Q, and showed that 
for the adequate description of hadron multiplicities the 
conservation laws should be modified, whereas for the 
description of hadron yield ratios the conservation laws 
are not necessary at all. 

Also here we analyzed the usual criteria for the chem- 
ical freeze-out and found that none of them is robust. 
Therefore we suggested a novel criterion of chemical 
freeze-out, a constant value of entropy per hadron num- 
ber equals to 7.18. Such a criterion is also supported 
by the different formulation of thermal model [2: and it 
evidences for the new physical effect which we called the 
adiabatic chemical hadron production. 

The performed analysis allowed us to find the restric- 
tions on the hard core radii that are imposed by the 
experimental data. Also we showed that although an in- 
clusion of the Lorentz contraction improves the fit qual- 
ity for any pair of baryon and meson hard core radii, but 
at the same time it has a small effect on the chemical 
freeze-out parameters and on K + /tt + ratio. 

In addition we introduced the surface tension in ther- 
mal model and made several global fits to find the pa- 
rameters of surface tension. Remarkably we found that 
the temperature of nil surface tension value depends on 
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the considered interval of the collision energy. Thus, if 
the low energy data are included into the fit, then the nil 
surface tension temperature is about 60 ± 5 MeV, while 
the data for the center of mass energies above 10 GeV 
lead to an essentially larger value of this temperature 
value T = 147 ± 7 MeV. The latter is a very intrigu- 
ing result since a very close estimate for the nil surface 
tension temperature was obtained recently within en- 
tirely different analysis of the quark gluon bag surface 
tension [35] . Therefore, it is possible that these two in- 
dependently obtained results, indeed, evidence that the 
(tri) critical temperature of the QCD phase diagram is 
between 140 and 154 MeV. 
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